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Understanding the interplay between spin and heat is a fundamental and intriguing subject. Here we report 
thermal spin injection and accumulation in CoFe/MgOAn-type Ge contacts with an asymmetry of tunnel 
spin polarization. Using local heating of electrodes by laser beam or electrical current, the thermally-induced 
spin accumulation is observed for both polarities of the temperature gradient across the tunnel contact. We 
observe that the magnitude of thermally injected spin signal scales linearly with the power of local heating of 
electrodes, and its sign is reversed as we invert the temperature gradient. A large Hanle magnetothermopower 
(HMTP) of about 7.0% and the Seebeck spin tunneling coefficient of larger than 0.74 meV K" 1 are obtained 
at room temperature. 

The tremendous power consumption and accompanying heat generation in current electronic devices 
requires alternative technologies to provide a solution for the energy issues and to realize the energy efficient 
electronics. Spintronics is a viable route for it, using the spin degree of freedom in addition to the conven- 
tional charge transport 1 . Recently another important ingredient, namely heat, appears on the stage. It has been 
reported that the temperature gradient and heat flow in ferromagnetic structures gives rise to a variety of spin 
related phenomena 2 . Understanding the interplay between heat and spin transport is a fundamental and intri- 
guing subject but also offers unique possibilities for emerging electronics based on the combination of thermo- 
electrics and spintronics. 

Especially in SC-based spintronics, the functional use of heat provides a new route to inject and control of spin 
in semiconductors (SCs) 3 . The generation of non- equilibrium spin populations (i.e. spin accumulation) in non- 
magnetic SC is a central issue of SC-based spintronics 1 ' 3 " 5 . Significant progress has been made on the spin 
accumulation in various SC systems by means of circularly polarized light 4 ' 6 , spin-polarized tunneling 7 " 13 , hot- 
electron spin filtering 1415 , spin-orbit interaction 1617 , or magnetization dynamics 1819 . 

Intriguingly, Le Breton et al. 20 have reported a rather different mechanism for the spin accumulation, in which 
temperature difference across a ferromagnet (FM)/oxide/SC tunnel contact can induce the spin accumulation 
(A/n) in SC via Seebeck spin tunneling (SST). It was found that the SST effect, involving thermal transfer of spin 
angular momentum from FM to SC without a tunneling charge current, is a purely interface- related phenomenon 
of the tunnel contact and governed by the energy dependence of its tunnel spin polarization (TSP) 20 . This provides 
a conceptually new mechanism for the generation of in SC as well as for the functional use of heat in spintronic 
devices 20 . The SST and thermal spin accumulation (Afi th ) inp-type SC (e.g.,p-type Si) have been intensely studied, 
using a Ni 80 Fe2o/Al 2 03/Si02/p-SC contact with the asymmetry in TSP of tunneling holes 20 . However, the 
counterpart of n-type SC still needs to be explored 21 . 

Here we report the thermal spin injection and accumulation in n-type Ge by employing CoFe/MgO/Ge 
contacts with the asymmetry in TSP of tunneling electrons. The Joule heating of Ge by electrical current or local 
heating of CoFe electrode by a laser beam gives rise to a temperature gradient across the tunnel contact and a 
consequent thermal spin injection into Ge. The thermally-induced spin accumulation was detected by means of 
the Hanle effect. The sign and magnitude of Hanle signals are analyzed using spin-dependent tunneling theory. 

Results 

Principle of the approach. Figures la and lb illustrate the thermal spin injection and accumulation process in a 
FM/oxide/SC tunnel contact for the case that the temperature of SC (T sc ) is higher than the temperature of FM 
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Figure 1 | Principle of the approach, (a) Schematic illustration of SST process in a FM/oxide/SC tunnel contact for the case that temperature of SC is 
larger than that of FM. (b) Spin-dependent density of states and its occupation for the tunnel contact with a hot SC and a cold FM. The representative 
profile of TSP vs. E for the CoFe/MgO/77-Ge contact is schematically illustrated on the middle. The (i) and (it) represent, respectively, the forward and 
reverse tunneling processes of electrons driven by the temperature difference, (c) Device geometry and measurement scheme. 



(Tfm) 20 - The temperature difference results in the unequal electron 
distributions in the FM and SC; the hot SC has a relatively larger 
number of filled (empty) states above (below) the Fermi-level (E F ) 
than those of the cold FM (see Fig. lb). Consequently, the electrons 
above E F mainly tunnel from the SC into FM (forward tunneling, 
indicated as (i) in Figs, la and lb), and the electrons below E F flow in 
the opposite direction (reverse tunneling, indicated as (ii) in Figs, la 
and lb). The important feature of the SST process is that the numbers 
of electrons tunneling in the opposite direction are the same, and 
there is no net charge current. Nevertheless, if the energy dependence 
of TSP are different for the forward and reverse tunneling (in other 
words, the Seebeck tunnel coefficient for majority and minority spin 
electrons are different), a spin accumulation (or Afi th ) can be induced 
in the SC 20 . Obviously, the sign of Afi th can be reversed when the 
temperate gradient is inverted; the cold SC has a relatively smaller 
number of filled (empty) states above (below) the Fermi-level (E P ) 
than those of the hot FM. 

Figure lc shows the device geometry and measurement scheme 
used in the present study. We fabricated a symmetric device consist- 
ing of three epitaxial CoFe(5 nm)/MgO(2 nm)//?-Ge tunnel contacts 
(a-c, 100 x 300 /nm 2 ) where the n-Ge channel is composed of a heav- 
ily P-doped surface layer ( rid ~ 10 19 cm _3 at 300 K) and a moderately 
Sb-doped substrate ( n^^ 10 18 cm _3 at 300 K). These contacts are 
separated by about 100 ^mfrom each other, which is much longer 
than the spin diffusion length. The magnetic easy axis of the CoFe 
contacts are along the [110] direction of Ge in parallel to the long axes 
of the contacts. Details of the structural and electrical characteriza- 
tions can be found elsewhere 22,23 . 

To effectively generate the temperature difference (AT=T Ge — 
T Co Fe) in the tunnel contact a, we used two different heating methods. 
For the SC heating (top panel of Fig. lc), we applied a heating current 
(heating) through the SC channel using two contacts b and c, which 
causes Joule heating and raises T Ge with respect to T CoFe (AT>0). 



For the FM heating (bottom panel of Fig. lc), the Au bond pad was 
heated using a laser beam with a wavelength of 532 nm and a max- 
imum power of 200 mW (note that the results obtained with the 
Joule heating of FM are shown in the Supplementary Information 
Note I). The laser beam spot (with a diameter of 5-10 |im and a skin 
depth of ~3 nm) on the 100 nm-thick Au pad is located at around 
300 |im from the one edge of the tunnel contact. The size of the Au 
pad is large enough to prevent the direct illumination of the Ge layer. 
A part of heat generated from the laser beam passes through the 
contact, resulting in T Co Fe > Tc e (AT < 0). 

In an open -circuit geometry, where the tunneling charge current 
(I tunnel) is zero > the measured voltage between the contact a and d is 
given by V = V t h + AVth 20 - The first term V tn is the thermo voltage 
maintaining zero net charge current (I tU nnei = Q) an d the second 
term A is the SST voltage due to the induced A/n th in the SC. 
The AVth can be detected by means of the Hanle effect 24 ' 25 . When 
we apply magnetic fields (B z ) transverse to the spins in the SC, 
the A/n th is suppressed via spin precession. This results in a voltage 
change (A V T h)> directly proportional to A/n th , with a Lorentzian line 
shape as a function of the B z . These two measurement schemes 
(Fig. lc) using the Hanle effect 24 ' 25 provide a concrete means of 
demonstrating the SST and resultant A/n th in the SC. 

Energy dependence of tunnel spin polarization in the CoFe/MgO/ 
n-Ge contact. Before conducting the thermal spin injection experi- 
ments, we have estimated the form of TSP as a function of £, which 
determine the sign and magnitude of induced Afi th , using the 
conventional three-terminal Hanle (TTH) measurement 7 ' 20 ' 26 . Instead 
of applying lh ea ting> a non-zero l tunne \ (or I ac ) is applied across the con- 
tact a and c of the CoFe/MgOM-Ge contact while the V is measured 
using the contact a and d in an applied magnetic field (see Fig. lc). 

Figure 2a shows the results obtained from the TTH measurements 
(up to 4 kOe) under perpendicular (B z , closed circles) and in-plane 
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Figure 2 | Energy dependence of tunnel spin polarization for the CoFe/MgO/n-Ge contact, (a) TTH measurements (up to 4 kOe) for the CoFe/MgO/n- 
Ge contact under perpendicular (B z , closed circles) and in-plane (B x , open circles) magnetic fields at 300 K. (b) Electrical Hanle signals (AVeh, normal 
and AVeh, inverted) as a function of an applied current (I ac ) at 300 K. (c) Corresponding spin RA products (AVeh, normal/ Jac and AVeh, inverted I lac) and (d) 
estimated TSP 2 for AV £H , normal I he with a bias voltage ( V e \ = V Co Fe - V Ge )> defined as V B = 0 - AV £H - 



(B x , open circles) magnetic fields at RT. Clear electrical Hanle signals 
(AVeh) with a Lorentzian line shape are detected at RT. It should be 
noted that the full spin accumulation (Afi) consists of the inverted 
Hanle signal (AVeh, inverted) in B x and the normal Hanle signal 
(AVeh, normal) in B 27 ' 28 . As shown in Figs. 2a and 2b, the electrical 
Hanle signals (AVeh) are significantly asymmetric with respect to the 
bias polarity, which is consistent with the previous work 23 . Both 
the AVeh, normal and AVeh, inverted increase linearly with increasing 
the reverse bias current (I ac <0, electron spin injection), but change 
slightly with increasing the forward bias current (I ac > 0, electron 
spin extraction). This asymmetry can also be seen in the middle panel 
of Fig. 2c, where the spin-RA products (AVeh /Jac) is plotted as a 
function of the electrical voltage (V e \ = V Co Fe — VGe)- It should be 
mentioned here that the obtained spin-RA value (2-3 kQ|im 2 ) is 
several orders of magnitude larger than the expected value from the 
existing spin injection and diffusion theory 5 . This discrepancy 
between experiment and theory in the TTH measurement has been 
consistently observed with many types of tunnel barrier and SC, as 
discussed in the Ref. 3. The origins of this discrepancy, other 
enhancement factors not yet incorporated in the existing theory, 
are still under investigation 29 . 

According to the spin injection and diffusion theory 5 , the 
AVeh, normal / Jac ^ proportional toy^y^^r^at a given tempera- 
ture (T), which depends on V d (note that, although the existing 
theory does not provide the quantitative agreement with the experi- 
ment, it provides the qualitative description of the TSP 2 variation 
with E). Here, the y^is the TSP corresponding to the detection of 
induced spin accumulation at the Ge interface, the y^ e is the TSP of 
the injected/extracted electrons, and the t s j is the spin lifetime. Using 
the AVeh, norma i /Jac vaiues (Fig- 2c) and effective (t 5 ^) values (not 
shown) extracted from the Lorentzian fit (black line in Fig. 2a), we 



plotted the TSP 2 (y^y,-/^ normalized) vs. V e \ in Fig. 2d. With the 
assumption of y d = y i / e , the variation of TSP (y) with E is then 
obtained asy(£)ocy 0 for£<£E (gray line in Fig. 2d) andy(£)ocy 0 
exp ( — (E — E F ) /0. 1) for E > E F (black line in Fig. 2d). The estimated 
profile of TSP vs. E is indeed similar to the schematically illustrated 
one on the middle panel of Fig. lb. This asymmetry in the TSP of 
CoFe/MgO/n-Ge contacts is essential to induce the large A/n th in Ge 
via the SST 20 . 

Detection of thermal spin accumulation in Ge. The SST and 

resultant A/n th in Ge were detected in the measurement geometry 
shown in the top panel of Fig. lc. While heating the Ge side 
(AT>0) as described above, the voltage change (AVth°c A/%,) i s 
measured as a function of the applied magnetic field (B z , B x ). 

As shown the A Vth~B z plots in Figs. 3a and 3b, significant normal 
Hanle signals (A Vth, normal) with the Lorentzian line shape similar to 
that of AVeh (see Fig- 2a) were observed at 300 K. The inverted 
Hanle signals (A Vth, inverted) 27,28 inB x (Figs. 3c and 3d), roughly half 
the magnitude of the A V T h, normah were also clearly measured. 

The magnitude of thermal Hanle signals (AVth) as a function of 
the heating ( U P to + 10 mA) is summarized in Figs. 3e-3h. In these 
figures, we can clearly see that the A Vth, normal and also A Vth, inverted 
scale quadratically with the I nea ting> and scale linearly with the 
heating power density (P heating)- The P heating is proportional to 
R Gell e ating/ L x eL v eL z e > where the R Ge is the resistance of the Ge heat- 
ing layer and tne VfL^ e L^ e is its volume. These results strongly 

support that the observed signals mainly come from the thermally 
driven spin accumulation, which scales linearly with the AT 20 ' 30 , in Ge. 

The sign of Afi th has been determined by a direct comparison 
with that of AVeh obtained from the TTH measurements for the 
same tunnel contact a 20 . The sign of AVeh, normal is negative for 
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Figure 3 | Detection of thermal spin accumulation in Ge with heating SC. Thermal Hanle signals ( A Vth, normal* A Vm, inverted) in applied magnetic fields 
(B z , closed circles; open circles) for heating currents (heating) of + 5 and + 10 mA at 300 K. (a)/(b) and (c)/(d) represent AVth, normal and A Vth, inverted 
for heating of + 5/ + 10 mA, respectively, (e) AVth, normal an d (f) AVth, inverted with the heating (up to + 10 mA), together with a quadratic fit. 
(g) AVth, normal and (h) AVth, inverted with theP heatingy together with a linear fit. 



negative V e \ (Vc 0 Fe — ^Ge < 0, electron spin injection) and positive for 
positive Vg/ (Vc 0 Fe — VGe >0, electron spin extraction) in Figs. 2(a) 
and (b). The AVth, norma/ in Figs. 3a and 3b has the same sign with 
the former AVeh, normal (V e i<0), which means that the Afi th pro- 
duced by the thermal spin injection has the same sign as 
the Afi el induced by the electrical spin injection. Given the positive 
TSP of epitaxial bcc FM/MgO(001) tunnel interfaces 31 ' 32 , the 
Afi th induced by the SST for T Ge > T Co Fe corresponds to the majority 
spin accumulation (Afi > 0) in the Ge. This is in agreement with the 
expected Afi th from the SST mechanism 20 : when the TSP of FM/ 
oxide/SC contacts (with positive sign) is constant below E F but decay 
above E F , as shown in Figs, lb and 2d, majority spins accumulate in 
the SC (Afi > 0) for T Ge > T Co Fe (see Figs, la and lb). We have calcu- 
lated that the induced A/n th is ( + )0.26 meV with the maximum 
Pheating (667 nW urn" 3 ) fromAn th = (-2e)AV TH /y 3,5 ' 20 , using the 
measured AVth, normal of (-)0.09 mV and the assumed TSP (y) of 



( + )0.7 for the epitaxial CoFe/MgO tunnel interface 31 ' 32 . This value 
should be considered as a lower limit for the Afi th , since we used the 
highest value of TSP. The SST theory 20 based on the free-electron 
model predicts the A/n th value of —10 |ieV for the AT of ~1 K. 
Considering the possible range of AT across the tunnel contact is 
from 0.15 mK to 350 mK (see the Supplementary Information Note 
II), the obtained signal is at least one order of magnitude larger than 
the predicted value, which was also observed in the Ni 8 oFe 2 o/Al 2 03/ 
Si0 2 /p-Si contact 20 . It is probably due to the limitation of the free- 
election model, the lack of detailed information about the energy 
dependence of TSP near the E F , and the ignorance of inelastic (mag- 
non-assisted) spin tunneling 20 . 

It should be mentioned that our system requires about 100 times 
larger Pheating^ 0 obtain a similar magnitude of AVth compared to the 
silicon-on-insulator (SOI) wafer based system 20 . This is mainly 
attributed to the different device structure. In our device structure, 
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the AT across the tunnel contact is small because a large part of the 
Joule heat produced in the heavily- doped Ge layer can flow away to 
the substrate. The calculated heat flow ratio (R^ nnel ) across the tunnel 
barrier in our device is more than one order of magnitude smaller 
than that in the SOI wafer based device 20 (see the Supplementary 
Information Note II). 

Sign reversal of thermal spin signal by reversing the temperature 
difference. Another important feature of SST and A^ is that the 
sign of the thermal spin signal is reversed when AT is reversed. To 
demonstrate this, we employed a laser beam for heating the FM 
instead of the Joule heating method. This approach allows us to 
heat the FM (AT<0) effectively and exclude the contribution of 
spurious effects such as the current-in-plane (CIP) tunneling and 
anisotropic magnetoresistance (AMR) on the Hanle signal (see the 
Supplementary Information Note I). 

Figures 4a and 4b show the obtained thermal spin signals as a 
function of perpendicular (B z , Fig. 4a) and in-plane (B x , Fig. 4b) 
magnetic fields with varying the laser power density (Pi aS er) at RT 
(AT<0). It should be noticed that we used the same contact a (see 
Fig. lc) in both cases of the SC and FM heating. Clear normal and 
inverted Hanle signals are observed and the amplitudes of both 
Hanle signals are gradually increased with increasing the P\ aser . As 
shown in the AVth -Piaser plots (Figs. 4c and 4d), the obtained Hanle 
signals (A Vth, normal and AVth, inverted) scale almost linearly with the 
Pia S er> indicating that the obtained signal is from the thermal spin 
injection and accumulation in Ge. 

We can see that the sign of A Vth is apparently reversed with 
respect to the SC heating case: the A Vth, normal (AVth, inverted) is pos- 
itive (negative) for T Co Fe > Tc e (see Figs. 4c and 4d) whereas it is 
negative (positive) for T Ge > T Co Fe (see Figs. 3g and 3h). This result 
clearly demonstrates another key feature of the SST and Afi th that the 
thermal spin signal reversed when AT is reversed. 

In addition, this sign of the Hanle signal allows us to exclude other 
possible origins for the spin signal, for example, due to the spin- 
polarized hot-electron injection and conventional Seebeck effect, in 



the laser-heating experiment. The thermal excitation by the laser 
heating can result in a current flow in the Au pad and FM layer. 
Nevertheless, this electrical current cannot be an origin of the obser- 
ved Hanle signal, since the sign of the observed A Vth, normal (posi- 
tive sign, see Fig. 4a) for the FM heating (T CoFe > T Ge ) by the laser 
beam is opposite to that AV^ nor mai (negative sign) expected by 
the injection of spin-polarized hot-electron current into the Ge 
(see the bottom panel of Fig. 2a). This strongly suggests that the 
observed Hanle signal does not come from the electrical current in 
the injector. 

Discussion 

For a quantitative analysis, we have calculated the Hanle magne- 
tothermopower (HMTP) (or Hanle magneto-Seebeck ratio 
(SHank)) 30 defined as the relative change of charge Seebeck coefficient 
(S) due to spin accumulation as follows: 



V\ 
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where G is the electrical tunnel conductance, L is the thermoelectric 
tunnel conductance, S 0 is the charge thermopower at zeroA^, 
andSssjis the SST coefficient. Figure 5a shows the calculated 
HMTP (orS H anie) as a function of the Pi aser . The HMTP value 
of AVth normal / Vth is ~ ( + )7.0% and remains almost constant with 
varying the P/ a5er . It is important to note that these values are 
—70 times larger than the maximum value (—0.1%) of Hanle mag- 
netoresistance (HMR), defined as AV^, normfl //y e /, obtained by the 
electrical spin injection using the same contact a (see Fig. 2). This 
experimental result supports the theoretical proposition 30 that the 
thermal spin injection is more efficient than the electrical spin 
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Figure 4 | Detection of thermal spin accumulation in Ge with heating FM. Obtained thermal spin signals under (a) perpendicular (B z ) and (b) in-plane 

(B x ) magnetic fields as a function of the laser power density (Pi aser ) at RT in the case of heating the FM (AT < 0). (c) A Vth, normal and (d) AVth, inverted as a 
function of the laser power density (P/ a5er ), together with a linear fit. 
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injection for the tunnel contacts with the strong asymmetry in TSP 
and the moderate thermal interface resistance. 

To estimate the value ofSsrr, defined asA/^/AT 30 , the infor- 
mation of AT across the CoFe/MgO/Ge tunnel contact is required. 
We carried out the combined analyses using the one-dimensional 
(1-D) heat flow model 20 and simulation using a commercial finite 
element software (COVERTOR) to determine AT (see the Supple- 
mentary Information Note II). Both the 1-D heat flow model and 
commercial finite element simulation provide a similar magnitude 
of AT. 

Figures 5b and 5c show the representative results of the simulated 
temperature distribution in a two-dimensional (2-D) cross -section 
and temperature line-scan cross the tunnel contact when heating the 
Ge with the maximum P heating of 667 nW |im" 3 . Here we used a 
reasonable value of 0.04 W m" 1 K" 1 (or a Gf h g ° eff of 2X 10 7 W 
m -2 y^~ 1 ) for the 2 nm-thick MgO barrier. The possible range of AT 
across the tunnel contact is from 0.15 mK to 350 mK, and a reas- 
onable AT is an order of 10 mK (see the Supplementary Information 
Note II for the detailed discussion). Using Afi th = 0.26 meV and AT 
= 10 mK, the Sstt value is estimated to be 26 meV K" 1 for the CoFe/ 
MgOM-Ge contact. If we use the largest A Tof 350 mK, the S S tt goes 
down to 0.74 meV K" 1 . These S S tt values are comparable to those of 
the Ni 80 Fe2o/Al 2 03/Si02/p-Si contact 20 but they are at least two 
orders of magnitude lager than the spin Seebeck coefficient ( — 3.8 
|aV K" 1 ) of the Co/Cu metallic system 33 , which may be attributed to 
the different spin injection/detection efficiency and the different 
measurement scheme. 

In conclusion, we have observed the thermal spin injection and 
accumulation in CoFe/MgOM-type Ge contacts using local heat- 
ing of electrodes by laser beam or electrical current. We dem- 
onstrate that the magnitude of thermally injected spin signal 
scales linearly with the power of local heating of electrodes, and 
its sign is reversed as we invert the temperature gradient. Based on 
the Seebeck spin tunneling theory, we have obtained a large Hanle 
magnetothermopower (HMTP) of about 7.0% and the Seebeck 
spin tunneling coefficient of larger than 0.74 meV K" 1 in the 
CoFe/MgOM-Ge contact at room temperature. This is attributed 
to the strong asymmetry in the TSP of the contact. Our experi- 
mental results demonstrate that the thermal spin injection is an 
efficient way of spin injection in SC. 



Methods 

Device fabrication. The single crystalline CoFe(5 nm)/MgO(2 nm)/rc-Ge(001) 
tunnel structures were prepared by molecular beam epitaxy (MBE) system with a base 
pressure better than 2X 10" 10 torr. Prior to the deposition of CoFe/MgO layers, the 
composite n-Ge consisting of a heavily P-doped surface layer ( rid ~ 10 19 cm~ 3 at 
300 K) and a moderately Sb-doped substrate ( rid ~ 10 18 cm~ 3 at 300 K) was formed 
by ion implantation technique. All layers were deposited by e-beam evaporation with 
a working pressure better than 2 X 10~ 9 torr. The MgO and CoFe layers were grown at 
125°C and RT, respectively, and then the samples were subsequently annealed in-situ 
for 30 min at 300°C below 2X 10" 9 torr to improve the surface morphology and 
crystalllinity. Finally, the samples were capped by a 2 nm thick Cr layer at RT to 
prevent oxidation of the sample. The symmetric device (Fig. lc) consisting of three 
tunnel contacts with lateral sizes of 100X300/100X300/100X300 una 2 was prepared 
by using micro-fabrication techniques (e.g., photo-lithography and Ar-ion beam 
etching) for the thermal spin accumulation experiments. For the electrical isolation at 
the sides of the tunnel contacts, about 120 nm-thick Ta 2 0 5 (115 nm)/Al 2 0 3 (2 nm) 
layers were grown by sputtering technique. For the contact pads, 110 nm-thick 
Au(100 nm)/Ti(10 nm) layers were deposited. 
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